Inside an evaporating two-dimensional charged black hole 
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We investigate the inner structure of an evaporating charged black hole, within the context of 
semiclassical dilaton gravity in two dimensions. The matter fields are charged, allowing the evapo- 
ration of both the mass and charge of the black hole. We find that the semiclassical effects cause the 
inner horizon to expand (by a finite factor), rather than to shrink to a point singularity. Although 
this expansion is a quantum phenomenon, the overall expansion factor is found to be independent 
of the magnitude of the quantum terms in the effective theory. 
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The internal structure of charged black holes (BHs) 
has been the subject of many investigations over sev- 
eral decades. The Reissner-Nordstrom solution ad- 
mits a perfectly regular inner horizon (IH) inside the 
BH, as opposed to the spacelike singularity of the un- 
charged Schwarzschild solution. Classical perturbations 
are known to develope instabilities on the IH, which lead 
to a null curvature singularity 1]. However, this singular- 
ity is found to be weak and non-destructive, in the sense 
that its tidal effect on a physical object is bounded and 
small . A similar situation was observed in a spinning 
BH 0. 

When semiclassical effects are taken into account, the 
black hole evaporates. In the case of a charged BH there 
are two distinct semiclassical phenomena: The electro- 
magnetic Schwinger effect, and the gravitational Hawk- 
ing effect. The evaporation of a charged BH involves 
both phenomena. An important question is, how does 
this evaporation affect the internal geometry. The goal 
of the present paper is to address this problem within a 
simplified two-dimensional (2d) model. 

Callan et al j3| developed a dilatonic 2d model of an 
uncharged BH. This model as well as similar ones were 
subsequently investigated by many authors 0. Simi- 
lar models of a charged BH (yet with uncharged mat- 
ter fields) were later developed 0. Nojiri and Oda 
then considered a dilatonic model with a charged 
quantum matter field. In their model, however, the 
dilaton-electromagnetic coupling took a non-standard 
form, which led to undesired causal structure of the 
charged BHs already at the classical level. This prob- 
lem was later fixed in Ref. ^9^. 

The model |9| consists of gravity in two dimensions, 
coupled to a dilaton 0, electromagnetic field Fuv, and 
N charged matter fields. These fields were taken to be 
massless chiral fermions, following Nojiri and Oda 0. 
The dilaton-electromagnetic coupling was taken as e"^"^ 
(rather than the coupling e^'^ used in Ref. ,8]). The 
analysis in focused primarily on the evolution of the 
various fields outside the evaporating BH. In particular 
the evaporation rate of the BH's mass and charged was 
determined. In this paper we shall use the same model 
to explore the evolution of geometry inside the BH, with 



special emphasize on the dynamics near the IH. 

Recently Frolov, Kristjansson, and Thorlacius (FKT) 
[lo| developed a similar model of 2d dilaton gravity cou- 
pled to a charged matter field. The main difference is 
that in Ref. [iflj the matter field is bosonized before 
semiclassical corrections are calculated. Also this field 
may be either massive or massless. In the massless case 
the field equations are closely related to those of Ref. |^ . 
FKT evolved the field equations (for a massless field) 
numerically in order to explore the dynamics inside the 
evaporating charged BH. They concluded that the com- 
bination of the Schwinger and Hawking effects turns the 
IH into a spacelike singularity where (f> diverges to -|-c». 
We employ here the usual translation of the 2d dilatonic 
terminology to that of a four-dimensional (4d) spherical 
model, in which the area coordinate is mapped to e^^"^. 
Then in this language, the inner horizon was found in 
[lo| to shrink to a point singularity. 

Previously several authors studied numerically the 
backreaction of Schwinger's pair-production effect on the 
internal structure of a charged BH, in either 4d ^3 o^' 
2d gravity (13|. These analyses, however, did not incor- 
porate the Hawking effect and hence the BH did not evap- 
orate, a fact which drastically changes the global struc- 
ture. Also Balbinot and Brady |Q| studied the backre- 
action to the Hawking effect inside a 2d charged BH, but 
the Schwinger effect was not included, and the BH was a 
non-evaporating one. To the best of our knowledge Ref. 
[ifij i is the only work published so far which analyzes the 
inner structure of an evaporating semiclassical charged 
BH. 

In this paper we explore the field equations of Ref. 
analytically, and reveal the internal structure of a semi- 
classical charged BH which evaporates its both mass and 
charge. We find a behavior rather different than that re- 
ported in On approaching the IH, the semiclassical 
effects cause the area coordinate to expand (to a finite 
value), rather than to shrink to a point singularity. This 
semiclassical expansion is linear in the Eddington coor- 
dinates. These analytical results were also confirmed by 
numerical simulations of our field equations. The origin 
of the discrepancy between FKT's analysis and our re- 
sults still needs be understood. In this paper we briefly 
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present our analysis and describe the main results. A 
more detailed description will be given in a separate pa- 
per 

The classical action and the semiclassical corrections 
are given in Ref. (it is the same as in Ref. 
apart from the modified dilaton-electromagnetic cou- 
pling). Our notation here follows that of Ref. 0, ex- 
cept that we choose our units such that the cosmolog- 
ical constant A = 1. We express the metric in double- 
null coordinates, ds^ = —e^'^dudv. (Both u and v are 
future-directed, as usual, and the lines of constant u 
or constant v are respectively parallel to past or future 
null infinity.) The electromagnetic field is expressed as 
Fuy = {g/y/Q)e^^'^^''^q, where g > denotes the electro- 
magnetic coupling constant. The system thus consists of 
three unknowns, ^, p, q. We then transform to new vari- 
ables R = e-^^ and S = 2{p - cj)). The resulting field 
equations 0| include a set of three hyperbolic equations 
for the three unknowns i?, S, q: 

R^uv = (qVR^ - l)e^ - Kp^uv , (1) 
S^uv = -2{qyR^)e'^ + Kp^uy/R, (2) 
q,uv = -Kg\q/R^)e'' , (3) 

where K = iV/24 (representing the magnitude of the 
semiclassical effects), and p = {S — lni?)/2. In addition 
there are two constraint equations, which both take the 
form 

R,ww Fi,wS^w '-^■ww — 5 (4) 

where hereafter w stands for either u or v, and 

Tnjw = K[{p,^^ - pI^) + {Kgy'^q^^ + z^{w)\ (5) 

represents semiclassical fluxes of energy in the u or w 
directions. The functions z„(u) and z^(u) are determined 
from the boundary conditions. 

The classical equations are obtained by formally set- 
ting K = Q. Then T^w vanishes, and q is set to a con- 
stant. The general solution then forms a two-parameter 
set, parametrized by the mass m and the charge q, to 
which we shall refer here as the two-dimensional Reissner- 
Nordstrom (2dRN) solution. This solution is most easily 
expressed in diagonal coordinates 0. For the present 
analysis, however, we shall need the classical solution in 
the double-null coordinates 

A key function which appears in the classical solution 
is H{R) = R — 2m + q^ /R. We are considering here the 
non-extreme case, m > q > 0. Then this function has 
two zeros, at R± = m± {vn? — q^Y^^. Correspondingly 
the charged BH admits two horizons: an event horizon 
(EH) ai R^ R+ and an IH at i? = i?_ . The latter also 
functions as a Cauchy horizon (CH) for initial hypersur- 
faces outside the BH. 

For the analysis below we shall need the classical solu- 
tion in between the two horizons, namely at i?_ < R < 



in the " Eddington-like" gauge in which both R and 
S depend on the single variable x = u -\- v. The solu- 
tion consists of the two functions R{x) and S{x)^ which 
satisfy 

- -H{R) , R^, - HiR) . (6) 
R{x) is formally expressed through its inverse function, 




[l/H{R')]dR' . 



(The lower integration limit was set here to m for con- 
creteness.) Note that H is negative at i?_ < R < 
and correspondingly x tends to —00 at the EH and to 
+00 at the IH. In a classical collapse scenario the EH 
(IH) is located at w — > —00 (v — s- -l-oo) and finite v {u). 
We also define the two parameters 

K± = ^{dH/dR)n=R^ = ±[(m2 - 9^)1/2/^^] . 

Note that in our notation > and k_ < 0, and that 
K± depend solely on q/m. 

The domain i?„ < R < is naturally divided into 
three regions: (i) The near-EH zone, namely the region 
of very negative x {—x » 1), (ii) The near-IH zone, 
X » 1, and (iii) The central zone, namely the region 
of moderate x values. In both the near-EH and near- 
IH zones e'^ is exponentially suppressed, like 6^*^+^ and 
g2K_2; respectively. 

Consider next the subset of classical solutions with 
fixed ratio q/m = a. It depends on a single param- 
eter m. One can easily verify that all solutions in this 
one-parameter family are related by a simple scaling rule: 

R{x; to) = mRo (x) , S{x; to) = In m + 5*0 (x) , 

where Ro{x) = R{x; m = 1) and Sq{x) = S{x; to = 1) are 
the functions corresponding to unit mass (and q — a). 

As long as the BH is macroscopic (to >> K), the semi- 
classical effects outside the BH are weak in a local sense; 
Namely, in the neighborhood of each point the metric is 
that of a classical 2dRN solution weakly perturbed by 
the semiclassical fiuxes. These weak fluxes accumulate 
over long time scales, however, and lead to slow decay 
of TO and q. We denote the slowly-drifting values of m 
and q, measured along the EH, by meh{v) and qeh{v) re- 
spectively. In Ref. jsl we calculated the rate of change of 
these quantities, and analyzed their long-term evolution. 
It was shown that for g < 1/2 the ratio q/m approaches 
a constant 

a=yT^/(l-g). (7) 

Furthermore, if q/m is initially a, then this ratio is pre- 
served during the evaporation. In this case the mass and 
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charge evaporate at constant rates (apart from correc- 
tions of order K/irieh Setting w = at the moment 
when (by extrapolation) rrieh vanishes, we may write 

mehiv) = -f3v , qeh{v) = -aPv , (8) 

where (3 — Kg^m/ R^. 

In what follows we shall consider the case g < 1/2. For 
simplicity we shall also assume that the BH's initial q/m 
is the preferred value (0 (but our main conclusions do 
not crucially depend on this assumption 15] ). 

We turn now to discuss the evolution inside the BH, 
assuming that the latter is macroscopic (throughout the 
rest of this paper by "macroscopic" we mean i?_ » K , 
which also implies m» K). In fact the above statement 
about the slow semiclassical evolution also holds inside 
the BH at i? > i?_ (except at the very neighborhood 
of as we discuss below). We can then express the 
solution for i?, S", g in a perturbative manner, where the 
leading order is the classical 2dRN solution with slowly- 
varying parameters m and q — am. We shall refer to 
this approximation as the adiabatic approximation. In 
the range i?_ < R < R+ which concerns us here, this 
expansion takes the form 

R{u,v) = ■meh[Ro{u + v)+eRi{u,v) + 0{e'^)] , (9) 
S{u,v) = \n{meh) + So{u + v) + eSi{u,v) + 0{e'^) , 
q{u,v) = meh[a + eqi{u,v) + 0{£^)] , 

where e = K/meh and, recall, rrieh = —[3v. The field 
equations 1114(1 are then automatically satisfied at the 
leading order e". At the next order one obtains a 
set of linear field equations for the three new functions 
Fi(u,?;), where hereafter F will stand for the three un- 
knowns i?, S, q. The initial data for Fi at the EH are 
discussed in Ref. [l5| . 

For the present analysis we do not really need the spe- 
cific solution for the perturbation fields Fi. All that we 
need is to make sure that these fields are well controlled 
inside the BH (and hence the perturbation term cx e is 
negligible compared to the leading classical term). A 
close look at the field equations for Fi reveals ^3 that 
indeed these fields are well controlled in the near-EH and 
central zones; However, in the near-IH zone (namely at 
large positive u -\- v) Ri and qi evolve linearly in u and 
V (see below) and get large values. This linear drift for- 
mally invalidates the above perturbative scheme at large 
positive u + v. Since in this paper we are primarily inter- 
ested in the behavior near the IH, we shall develope now 
another perturbative scheme which will be valid at large 
positive u + V despite the linear drift in R and q. 

We first bring Eqs. 1(113(1 to a standard hyperbolic form 
by isolating R^uv and S^uv We obtain a system of the 
form 



with Aq = -Kg^qjR? and = 0. For F = R, S the co- 
efficients Ap, Bp are certain functions of their arguments 
(and K) which we write explicitly in Ref. For the 

present analysis we shall not need the _B-functions, and 
also it will be sufficient to express Api and As at the 
leading order K'^: 

Ar ~ q^/R^ - 1 , As^ -2gVi?3 . 

In Ref. ^3 we analyze the effect of the terms propor- 
tional to R^uR.v on the near-IH asymptotic behavior and 
find it to be negligible (note that vanishes at the IH 
at the leading order). Omitting these terms, the evolu- 
tion equations take the simple form 

F„„ = ^F(i?,'?)e^. (11) 

Our goal is to explore the asymptotic behavior of this 
system at large u + v. To this end we first consider the 
classical case (namely K = and q — const). From 
Eq. © it follows that S decreases at large x linearly, 
S w 2k_(m -I- v). Consequently e"^ becomes negligible, 
and Eq. 1(11(1 reduces to 

F.uv=0- (12) 

The above asymptotic form of S obviously satisfies this 
equation. The equation for R implies R = Ru{u) + 
Rviv)+const, and it is indeed satisfied with i?„ — Ry ~ 
and const = i?_ . 

Turn now again to the system 1(11(1 in the semiclassical 
case. The most crucial ingredient in the above classi- 
cal asymptotic behavior is the global decaying factor e"^. 
There is no reason to expect that this behavior will be 
sensitive to small modifications in the coefficients Ap. 
Thus, the system ((11(1 seems to admit a rather generic 
class of asymptotic solutions at large u + v, in which S 
decreases linearly (or alike) in u and v, leading to Eq. 
((12(1 . Then the asymptotic solution should take the sim- 
ple form 

F = Fu{u)+F4v) (13) 

(up to exponentially-small corrections). The claim that 
our system of evaporating BH indeed falls into this class 
of asymptotic solutions will be justified below, by deter- 
mining the actual functions Fu and Fy. 

The functions Fu,Fy must satisfy certain consistency 
conditions in order to enable this simple asymptotic be- 
havior. First, Su and must decrease sufficiently fast 
in u and v respectively (a linear decrease, as in the clas- 
sical solutions, is sufficient). Second, since the functions 
Ap diverge at small R (in fact Ap has poles at i? = 
and R — K), the functions Ru and Ry must have a form 
which guarantees that R will increase (or at least not 
decrease) with u and v. 

In fact, the analysis below indicates that the asymp- 
totic solution p3(l takes the more specific form 



F,™ = Ap{R, q)e^ + KBp{R)R,uR,y , (10) 



F ^ apu -\- bpv + cp , 



(14) 
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where ap,bp,CF are certain constants. In terms of these 
parameters, the above consistency conditions reduce to 
(i) as, bs < 0, and (ii) a^, bu > 0. The sign of Ug and bg 
is unimportant. Note that with such a hnear behavior, at 
large u + v the source terms in Eqs. l(TT|) are aU bounded 
by (u + w)^e*^°^"+''^"\ which indeed has a neghgible effect 
and justifies the approximation H12I13|I . 

The determination of the required functions i^^(w) is 
done by matching the near-IH approximation H12(l to the 
adiabatic expansion 0. This matching is done in the 
intermediate zone 

1 «u + v « R.^/K . (15) 

The inequahty u + v >> 1 guarantees that e'^ oc 
^2k-(u+v) I Qj^^ ^}jg near-IH approximation ifT^ is 
vaUd. The inequality u + v<< i?_ / K means that in the 
adiabatic expansion the terms Fi (which are cx K 
but undergo linear drift at large u + v) are unimportant, 
and the adiabatic approximation holds. A consistency 
condition for this matching is that in the intermediate 
zone F^u must be independent of v, and vice versa. If in 
addition F^^ is independent of w, then the linear form 
(fTIjl holds as well. 

The parameters as and bs are nonvanishing already 
at the classical level {as — bs — 2k_). The semiclassical 
effects lead to 0{e) corrections to these parameters, but 
these small corrections are not needed here. On the other 
hand, a^j, 6^^, a^, 6g all vanish at the classical level, and 
our goal is to calculate these four parameters at order e^. 
It will be sufficient to calculate cp a,t order £°, because 
the linear drift in F is not sensitive to these constants. 

Consider first the semiclassical evolution of q. Inte- 
grating the field equation l(TT|l for F = q with respect to 
u or V, yields 

/oo 
{qe'^/R') dw , 
-oo 

where w is the null coordinate other than w (e.g. w = v 
when w stands for u) , and the integration is done along a 
line of constant w. Since we only need g^^, at order K^, it 
is sufficient to evaluate the integrand at the leading order 
namely it can be taken as ae^° / Rq. Taking the upper 
integration limit as oo is justified because of the exponen- 
tial decay of e'^ at large u + v. The integration is then 
straightforward and the integral becomes q/R- — q/R+- 
[Note that unlike q and R±, q/R± is independent of rrieh 
and only depends on a.] Noting that q'^u'^^ = (e.g. from 
regularity) and q'^v'^^ = qeh.v is —a/? = —Kg^q/R^, we 
obtain at the intermediate zone 

q^u = Kc/{q/R+ - q/R-) = a,, = -Kg^qjR^ = b^. 

Both quantities are independent of v and u, confirming 
the linear form (|14|l for q. 

Next consider the semiclassical evolution of R. Here 
again we evaluate R,u and R^^ at the intermediate zone. 



To this end we use the constraint equation Q). Since we 
only need R^w at order ^ it will be sufficient to cal- 
culate TujM) at its leading order — namely that of a test 
charged field in a classical 2dRN background. This test- 
field T^u, is to be evaluated at the near-IH zone of this 
2dRN geometry, where p^ww — and p^w — k-. Also the 
boundary conditions imply z„ = and z^, = 0. One 
then finds that T^u and Ty^ are just constants. Equa- 
tion Q then yields = T^wl S,w (the corresponding 
homogeneous solution R ^] = e^ may be ignored at large 
u + v). Substituting = 2k_ we find 

R^u ^ K[g\q/R- - q/R+)^ + (4 - n^_)] /2n_ = a„, 

and 

Rv = K[igq/R^)^ - kI] /2k_ = bn . 

Again, both and R^v are independent of v and u, 
which justifies the linear form (|14|) . 

As was discussed above, our near-IH approximation 
will only be valid if both an and 6^ are positive. Noting 
that Eq. Q implies = g, one obtains 

an =. Kg\R+/R^ + R-/R+) , b^ = Kg\R+/R_) . 

Both parameters are positive as desired. 

The constants cb, and Cq are calculated from Eq. H14|l 
by setting R ~ i?!'' and q ~ qeh at the intermediate 
zone, where denotes the value R± corresponding to 
iJiehiv). One obtains 01 Cq = c^^ = 0. 

The global picture of the region between the two hori- 
zons as emerges from our analysis is fairly simple. Since 
we have set u = at the end-point of evaporation, the 
CH of the evaporating BH is the line w = 0. The EH is 
located at it = — oo, just as in the non-evaporating case. 
V is negative and increasing along the EH, and u is posi- 
tive and increasing along the CH (starting from u « 0). 
In both the near-EH and the central zones the spacetime 
is well described by the adiabatic solution @, in which 
the mass and charge slowly evolve with v. In particular, 
R (like q) decreases linearly along the EH. In the near- 
IH zone all dynamical variables drift linearly in u and 
V. Along the CH R (like —q) grows linearly, R ~ a^w, 
starting from i? ~ at u ~ and subsequently growing 
to macroscopic values. 

In considering the near-IH behavior of q, one recognizes 
that both aq and bq are negative. In fact the charge flips 
its sign before the CH is approached, and along the CH 
q ~ aqU < 0. 

The positivity of aR and bn means that the effect of 
the semiclassical terms near the IH is to expand the latter 
(namely to increase R) , rather than to shrink it to a point 
singularity i? = 0. Consider an infalling observer which 
crosses the EH at some v = vq in the macroscopic domain 
(namely —vq » 1) and takes a typical orbit which in the 
locally-2dRN background would hit the IH at some finite 
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u. The observer approaches the central zone at v vq 
and M « Wo = —vq. Subsequently on heading towards 
the IH V increases remarkably but u is essentially frozen 
at M « mq. In the first part of the journey, during the 
adiabatic phase, R steadily decreases from i?+ to i?_ 
(corresponding to the entrance mass mph{v{))) just as in 
a 2dRN spacetime. But subsequently, on approaching 
the very neighborhood of the CH, R will start increase 
according to 

R ~ aR\vQ\ + hRV ~ Rt^ivo) + hR{v - vq) . 

Note that this phase of increasing R lasts an enormously 
short proper time, due to the huge blue-shift factor near 
the IH. The final value at the CH is i? ~ ai^jwol- A 
simple calculation shows that this value is even greater 
than Rf{vo), by the factor + > 1. 

Note that this final R value at the CH is indepen- 
dent of K; Namely, although the near-IH increase in R 
is a quantum phenomenon, its magnitude only depends 
on the macroscopic parameters irieh and qeh , not on the 
amplitude of the quantum terms in the effective action. 
This has a simple explanation: The rate of the near-IH 
semiclassical increase in i? is oc K, but at the same time, 
since P (x K, for a given entrance mass mehivo) the evap- 
oration time (from meh(vo) to complete evaporation) is 
oc 1/K, hence their product is independent of K. The 
same considerations apply to the semiclassical drift in q, 
and one finds that the final charge at the CH is 2k_ times 
the entrance value Qehivo), independent of K. 

We also performed numerical simulations of the full 
nonlinear hyperbolic system [using the form HlOfl 

with the exact functions Ap,Bf]. We started from ini- 
tial data corresponding to Eq. ||SJ) along the EH. These 
simulations confirm the above analysis. In particular, all 
quantities F where found to evolve linearly in u and v 
in the near-IH zone, with R increasing and q decreas- 
ing. The numerically computed quantities F„ and F^y 
all agree with the above analytically calculated parame- 
ters ap and bp to within a few percents. More details of 
the numerical simulations will be presented in Ref. [isj . 

The rapid change in i? as a function of the observer's 
proper time in the near-IH zone rases the question of 
whether the semiclassical theory is valid in this region. 
Scalars made of gradients, such as e~'^''(j)^u4',v, get enor- 
mously large (though finite) values on approaching the 
CH. Yet the effective gravitational coupling constant e* 
remains small (except near u — v — 0). Nevertheless 
the present analysis makes one point very clear: When 
the semiclassical theory comes to its limit of predictabil- 
ity (either at v = or at an earlier stage where gradient- 
based scalars become large), R^v is positive rather than 
negative. This fact may have important implications to 
the issue of extension of geometry beyond the CH, which 



we shall not address here. 

Finally we comment on the relevance of our results to 
the real 4d world. A key ingredient in the above 2d anal- 
ysis is the negative sign of Tuu and Tyy of a test field 
near the IH. In the case of a 4d spherical charged BH it 
is possible to bring the semiclassical field equations to a 
form very similar to the above 2d model. This similar- 
ity strongly suggests that in the 4d case too, the crucial 
factor which will determine whether the CH expands or 
contracts will be the sign of T„„ and Tyy near the IH, 
for a test quantum field in the Reissner-Nordstrom ge- 
ometry. Unfortunately in the 4d case there is no simple 
expression for T^i,. Nevertheless certain methods were 
developed for calculating T,,,^ in the Reissner-Nordstrom 
geometry outside the EH Il8]. The present analysis em- 
phasizes the importance of extending the calculation of 
Tpi, to the interior of the 4d Reissner-Nordstrom BH, 
particularly to the neighborhood of the IH. 

I am grateful to Larus Thorlacius for warm hospitality 
and fruitful discussions. I also thank Joshua Fcinberg 
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